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In a recent publication [B. Schiitte, S. Bauch et al., accepted for publication in Phys. Rev. Lett. 
2012] it was demonstrated both, experimentally and theoretically, that Auger electrons are subject 
to an energetic chirp if a three-body-interaction dubbed post-collision interaction (PCI) is involved. 
Here, we extend previous theoretical work and give a detailed analysis of field-assisted PCI based on 
numerical solutions of the time-dependent Schrodinger equation, extensive Monte-Carlo averaged 
molecular dynamics simulations and analytical theory. The dependence on various streaking and 
excitation conditions is investigated, and we discuss, how these findings may help to improve XUV 
pulse characterization as well as understanding of ultrafast atomic processes. 

PACS numbers: 32.80.Aa,32.80.Hd,32.30-r.78.47.J- 



I. INTRODUCTION 

The progress in the creation of phase stabilized laser 
systems and the generation of short and ultrashort pulses 
in the ultraviolet (UV) and extreme ultraviolet (XUV) 
regime [I] |2] allows nowadays for the observation of 
electronic processes on the femtosecond and even sub- 
femtosecond timescale in a time-resolved fashion [3l |4] . 
Fundamental investigations include the mapping of the 
oscillating electrical field of a laser ^, electron tunnel- 
ing in strong fields [6] or the direct observation of Auger 
decay in the time domain T]. Recently, processes down 
to a duration of several tens of attoseconds have been 
demonstrated to be resolvable [5]. 

The major tool for the observation of fast processes 
since early days in physics is the streak- camera. Fol- 
lowing its first mechanical realization by Wheatstone 
1834 [5] with /xs resolution, nowadays the sub-picosecond 
regime can be accessed with classical optoelectronic se- 
tups [ini [H] . To overcome the mechanical and electronic 
barriers for switching times, the answer was found in elec- 
trodynamics, leading to a setup called light-driven streak 
camera [12l HI]: Here, the temporal deflection of elec- 
trons is realized by the time varying vector potential of 
a laser field and the triggering of the process is done by 
ultrashort ionization through attosecond XUV pulses in 
pump-probe setups. A possibility to reach the zeptosec- 
ond regime in ultrahigh fields has recently been proposed 
[Hj theoretically. 

An important application of the light-field driven 
streak camera is the characterization of (X)UV pump 
pulses in the femtosecond 15, and sub-fs regime P3]- By 
means of photoionization of rare-gas target atoms the 
XUV pulse properties, such as duration, substructure 
and chirp, are imprinted on a photoelectron distribution. 
A time- varying streaking field deflects these electrons and 
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maps the temporal properties to a measurable energy 
spectrum. Therefore, this procedure strongly relies on 
the precise knowledge of the photon-to-electron conver- 
sion and, with that, a method to extract the temporal 
pulse properties from the streaked kinetic energy spec- 
tra of the electrons. While for solely photoelectrons this 
mapping is agreed to be understood for atoms [131 US] 
and atoms on surfaces [121 HT] , the situation strongly dif- 
fers if Auger decay is involved. 

The radiationless decay of resonances, flrst described 
by Lise Meitner 1922 [TH] and Pierre Auger 1925 [19], 
is a fundamental correlation-driven many-body effect in 
quantum mechanics, covering atoms, solids, quantum 
dots and molecules. The photoexcited inner-shell-hole 
is subject to spontaneous decay, transferring its energy 
to an outer-shell electron, that leaves the ion with its ex- 
cess energy. The result is a doubly charged ion and two 
correlated electrons in the continuum. Due to the spon- 
taneous character of the hole decay, the Auger electron 
cannot carry information about the pump pulse. Nev- 
ertheless, in pOj it was demonstrated, utilizing the THz 
streak camera setup, that the energy of the Auger elec- 
tron depends on its release time, i.e. it carries an ener- 
getic chirp, if the Auger electron is faster than the preced- 
ing photoelectron. This could be verifled using two inde- 
pendent experiments involving XUV photons from (i) the 
free electron Laser FLASH at DESY (Hamburg) and (ii) 
a higher-harmonics generation (HHG) source [2T]. While 
for (i) the ionizing pulse has a complicated structure, 
both in time and energy, for (ii) chirp-free pulses with 
rather well-deflned properties are expected. Neverthe- 
less, it was established experimentally that for both, (i) 
and (ii), the Auger electron's chirp is present and has 
qualitatively the same properties. 

The authors of [20] identified post-collision interac- 
tion (PCI) as the responsible mechanism for the observed 
chirp, utilizing extensive molecular dynamics (MD) sim- 
ulations as well as an analytically solvable model. PCI 
is a process, where a fast Auger electron can catch up 
with the slower photoelectron, which leads to a dras- 



tic change of the screening of the ion's charge. This 
manifests itself in an energy exchange: the photoelec- 
tron loses energy (increased binding), whereas the Auger 
electron is correspondingly accelerated (the binding po- 
tential becomes shielded). Obviously, the net amount of 
transferred energy depends on the distance from the ion, 
the closer the overtaking happens the stronger the effect. 
Although widely discussed in the literature [52H17] the 
consequences of PCI for the temporal energy distribu- 
tion remained unexplored. In this paper, we extend the 
theory presented in [5D| and give a detailed description of 
field-assisted PCI (FA-PCI) using quantum and classical 
simulations as well as analytical theory. 

The paper is organized as follows: In the first part. 
Sec. Inj we demonstrate the presence of a chirp on the 
Auger electron energy by solving the time-dependent 
Schrodinger equation (TDSE) for model systems and 
support the idea of PCI being the responsible mecha- 
nism. To overcome the model character necessary for 
quantum calculations, we develop in Sec. |III| a classical 
simulation technique based on Monte-Carlo (MC) aver- 
aged molecular dynamics (MD). In Sec. IV extending 
the model of |2D], we present an analytical theory for 
the Auger line shape in the presence of a slowly vary- 
ing streaking field including PCI effects and compare to 
TDSE as well as MC-MD simulations. In Sec. |V| we in- 
vestigate the influence of various pulse parameters and 



show, how the measurement of the PCI-induced chirp 
may help to improve the pulse characterization capabili- 
ties of light-field-driven streak camera setups. The paper 
closes with a comment on recent experiments and an out- 



look on future investigations in Sec. VI 



II. 



QUANTUM THEORY OF LASER-ASSISTED 
AUGER DECAY 



Starting point for the description of laser-assisted 
Auger decay (LAAD) on a quantum mechanical level is 
the time-dependent Schrodinger equation (TDSE). How- 
ever, full quantum calculations of autoionization involve 
two or more electrons, which limits them to model stud- 
ies, see e.g. |5S], or Helium ^29, on very short time scales. 
To overcome this "brute-force" approach we use a gen- 
eralization of Fano's theory [30] to the time- dependent 
case developed in ^31n55] and references therein. The 
notations follow [33], an analogous derivation based on 
quantum field theory can be found in I34j . A similar 
theoretical approach to time-resolved Fano resonances is 
developed in [351 ISS] ■ 

The time evolution of the outgoing photoelectron after 
excitation with an XUV pulse is governed by a set of 
coupled TDSEs (throughout atomic units, TTig = |e| = 
h — Aireo = 1, are used): 



^^^Mr,t) = (n^ir) -EA + ^[kA- ^^(i)]^ - z^i(i)') 0.(r,i) + VMr,t) 



(1) 
(2) 



Eq. (fTl) describes the photoelectron excited from the ini- 
tial orbital (f)Q with energy eg by a laser pulse Exit). It 
moves in a potential of a singly charged ion, included in 
Hi{r), and the streaking field ELit). In other words, 
0d(i) describes the photoelectron before decay of the res- 
onance with decay constant T a- After Auger decay with 
excess energy Ea, being the energy difference between 
the outer shell electron and the core hole, the photoelec- 
tron's movement in the potential of a doubly charged 
ion, contained in H2^ is governed by Eq. (l2|) coupled to 
Eq. (IT]) via the Auger decay matrix element V ^ which 
is assumed to be constant in energy and space |33| . By 
setting Hi = H2, post-collision effects due to changed 
screening of the ion's charge can be artificially excluded 
from the calculations. 



In fact, Eq. (pi) represents a set of equations for all 
possible energies of the Auger electron e = k\/2. The 
vector potential Al (t) associated with the electrical field 



ELit), 



AUt) = - f dTE^iT) , 

J —OG 



(3) 



is chosen to vanish for long times. 

Both laser pulses are linearly polarized in z direction 
with Gaussian envelopes and coupled to Eqs. ([I]) and ^ 
in dipole approximation. The streaking pulse with dura- 
tion tl, phase shift (p^ and frequency uj^ is centered at 
zero, 

EL{t)=e,Elexp( -^j cos [LULt + if l] ■ (4) 

The XUV pulse is delayed by tx with photon energy ojx 
and duration tx, 

Exit) = e,E°x exp (- ^^^^^ ) cos [oJxit - tx)] ■ 

(5) 
Note: throughout this paper, all pulse durations are 
given as full width at half maximum (FWHM) and will 



be denoted by tx.l = 2\/21n2fjc,L- The model, Eqs. ([!]) 
and ([2]), has been successfully applied to the recapture of 
photoelectrons due to PCI [371 [3H] and to (angle-resolved) 
sideband structures in LA AD [221110], which appear if the 
duration of the pump pulse is comparable to or longer 
than the period of the streaking field. In this work, we 
use Tx ^ l/'i^L required for streak cameras. 



A. Simplifications 

Up to now, the above-mentioned previous works con- 
sidered short pulses in the (sub-) fs regime involving 
infrared (IR) streaking pulses. The characterization of 
pump pulses longer than 20 fs, as they are produced e.g. 
by free electron lasers, requires deflecting fields based on 
THz radiation [TSl [20] and, therefore, requires the prop- 
agation of Eqs. (IT]) and (l2| over a duration of several pi- 
coseconds. In order to describe the involved processes on 
a time-dependent quantum mechanical level drastic sim- 
plifications are needed to keep the computational costs 
manageable. 

As a first step, we restrict our investigations to a one- 
dimensional (ID) version, i.e. consider wave functions of 
the form (f>d{x, t) and 4>e{x, t), neglecting any angular mo- 
menta and distributions. This leads to the model Hamil- 
tonians Hi (x) and H2 {x) which are chosen to account for 
the correct asymptotics of the binding potentials of the 
remaining ion. 



H,{x) 






Z, 



Vx2-t-, 



(6) 



with Zi — —1 and Z2 — —2. The Coulomb singularity 
appearing in ID systems has been regularized in a stan- 
dard procedure by k, e.g. [28l HH [42], assuring a finite 
binding potential at the position of the ion. 

Still, to keep track of the photoelectrons traveling with 
25 to 80 eV in the continuum, enormous computational 
grids are needed. To overcome this point, we introduce a 
scaling procedure of all relevant temporal quantities by 
a factor 7, which maps the (not-manageable) physical 
system to a smaller-sized analog, which can be tackled 
by the quantum simulations: 



^*A = I'^A, ^l = l^L, T*x = 



Tx 

1 ' 



TL_ 

1 



(7) 



In order to keep the relevant streaking conditions compa- 
rable, the intensity of the streaking field is chosen such 
that the ponderomotive potential Up = Eq/AloJ^ of the 
streaking field is kept constant when 7 is varied. The 
influence of this scaling procedure is discussed in detail 
below. 



B. Solution of Eqs. ([TJ and ([2J| 

We solve the ID analogs of Eqs. (II]) and ^ employ- 
ing a finite-element discrete variable representation (FE- 
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FIG. 1. (color online) Auger line shapes obtained by solving 
Eqs. ([l| and ([2| for a series of time delays tx- Shown is the 
case including PCI (blue solid lines, blue area) and the case 
neglecting PCI (red dashed lines), i.e. H2 = -ffi in Eqs. ([l] 
[2| for a 16.9 fs single-cycle (ipL = 7r/2) streaking pulse with 
a frequency of lol = 33 THz and Up = 98.8 meV. The 2.83 fs 
XUV pulse has a photon energy of 91 eV. The Auger de- 
cay constant was set to Fa ~ 950 meV at an Auger energy 
of Ea = 34.27 eV thus resembling THz-streaking of the Xe 
NOO transition scaled by a factor of 7 = 10 (see text for 
details). Note: the lines neglecting PCI have been shifted 
towards higher energy by 0.5 eV for better comparison. 



DVR) [331 m] and an independent finite-difference based 
method on large spatial grids allowing for the propaga- 
tion of several tens of fs without reflections at the grid 
edges. All considered observables have been carefully 
checked for convergence with respect to the numerical 
discretizations. 

Throughout this paper, two transitions motivated by 
the experiment are considered, the NOO transition in 
Xenon and the MNN transition in Krypton ^20,. Let us 
start with Xenon. For that, the eigenstate of Hi (k — 
0.1935) with a ground state energy of Ep — — 66eV is 
used for XUV excitation with a photon energy of wx — 
91 eV which corresponds to a kinetic energy of wx+Ep — 
25 eV for the photoelectron. The Auger electron energy 
is chosen to match Ea = 34 eV, being faster than the 
photoelectron wave packet and thus giving rise to PCI 
effects. An example of the resulting Auger electron line 
shapes for a full scan of time delays tx is shown in Fig. [T] 
for a streaking with 33 THz and Up — 98.8 meV and a 
pump pulse duration of 2.83 fs, thus scaled by a factor 
of 7 = 10 in comparison with the THz streak camera 
in [15. Analogously, the atomic parameters, matching 
the Xe NOO transition, are scaled by the same factor 
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FIG. 2. (color online) FWHM of Auger lines (right axis) ob- 
tained by solving Eqs. ^ and H for different scaling factors 
7, cf. Eq. ^, of the Xe NOO transition in a 3.3 THz streak- 
ing field keeping Up = 98.8 meV constant. The natural line 
width Ta has been subtracted for each set of parameters for 
better comparison. The maximum of the line (proportional to 
—A{t), gray line with filled circles labeled by -Emax, left axis) 
and the case neglecting PCI (black dashed line) are shown for 
7 = 10. For 7 = 10, parameters are the same as in Fig. [l] 



according to Eq. ([7|. 

Each individual Auger line is shown for two cases: in- 
cluding PCI (blue solid lines) and neglecting PCI (red 
dashed lines). For both cases, the typical streaking pic- 
ture of the time-dependent momentum transfer arises, 
with a general shift of the PCI result towards higher en- 
ergies (0.5 eV), which is compensated in Fig. [Tlfor better 
visibility. Careful inspection of the line shapes reveals, 
that for falling slope of A{t) {tx < 0) the lines are higher 
and of smaller width than for the case of rising slope of 
A{t) {tx > 0), cf. pci vs. nopci curves. Note, that the 
energy shift is proportional to —A{t). The lines corre- 
sponding to the case without PCI have the same height 
and width for positive and negative time delays. 

This observation already indicates a chirp in Auger 
electron emission, i.e. a time-dependent variation of the 
energy of the Auger electron manifesting itself in an 
asymmetry with respect to the direction of the slope of A. 
In the following, this result will be investigated in detail 
and the underlying physical mechanism will be identified. 



C. Analysis of the TDSE results 

Let us first discuss the influence of the scaling proce- 
dure (ItI), shown in Fig. [2J Wc point out, that each value 
of 7 corresponds to a certain physical system, but our aim 
is to describe experiments based on the Xe NOO transi- 
tion. The width displayed in Fig.[2]is extracted from line 
shape data by interpolation utilizing cubic splines and 
subsequent finding of the maximum and the correspond- 
ing FWHM. For better comparison, the a:-axis is shifted 
by the Auger decay time T^^ for each data set and the 
width was modified by -^/cr^ — r\ to account for the dif- 
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FIG. 3. (color online) FWHM of Auger lines (right axis) 
for different photon energies ujx of the XUV pulse of 28 fs 
duration. All other parameters are the same as in Fig. [2] for 
the case 7 — 10. The energy of the maximum of the Auger 
line and the case neglecting PCI (black dashed line) are given 
for ljjx = 91 eV. The time dependence of the electrical field 
E{t) is sketched by the gray dotted line. 



ferent natural line width in each set of parameters. The 
first observation is a strong asymmetry in the FWHM for 
all values of 7 with respect to the slope of A{t). Note: 
the displayed curve labeled i^^max shows the energy cor- 
responding to the maximum of the Auger line, which is 
proportional to —A{t). Approaching the physical sys- 
tem of Xe NOO (7 = 1), the asymmetry gets smaller, 
but is still present for the smallest considered value of 7. 
For comparison, also the case neglecting PCI is shown 
for 7 = 10, where no such asymmetry is observed and 
the typical chirp-free streaking behavior |13j is retrieved: 
largest width (and corresponding time resolution of the 
streak camera) occurs at maximum slope of A{t). Note: 
for single-cycle pulses used here, this does not coincide 
with zero transitions of A{t) (electrical field maxima). 
At the maximum of |^(0|, as expected, a pronounced 
minimum can be observed. 

Since the PCI effect originates from the changed 
screening of the remaining ions' charge during overtaking 
of the photoelectron by the Auger electron, it strongly de- 



pends on the velocity of the photoelectron, cf. Sec. IV C 



Therefore, the observed asymmetry should be more pro- 
nounced for slow photoelectrons, where the overtaking 
happens in close vicinity to the ion |20j , and should van- 
ish for fast ones, where the Auger electron cannot catch 
up with the photoelectron. Fig. [3] shows the FWHM 
of the Auger line of the Xe NOO transition (7 = 10) 
for a set of photon energies ujx- As is clearly seen, the 
strongest asymmetry is observed for slow photoelectrons 
(green curve with triangles) whereas the increase of the 
photoelectron's energy leads to a decrease of the observed 
asymmetry in the FWHM and approaches the case ne- 
glecting PCI (black dashed line) , thus supporting the idea 
that PCI is responsible for the energetic chirp in Auger 
emission. We note that, although for ujx — 126 eV rather 
fast photoelectrons (47 eV in comparison to 35 eV Auger 
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FIG. 4. (color online) FWHM of Auger line (right axis) for dif- 
ferent XUV pulse durations Tx^ ■ The maximum of the Auger 
line (proportional to ~A{t), filled circles) is shown for 28 fs. 
The time dependence of the electrical field E{t) is sketched 
by the gray dotted line. Parameters are the same as in Fig. [2] 
for 7 = 10. 
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FIG. 5. (color online) The same as Fig. |4] but for different 
ponderomotive potentials Up of the streaking field at a fixed 
XUV pulse duration of 28 fs (7 = 10). The graph of the 
maximum of the line corresponds to Up = 98.8 meV. 



electron energy) are emitted (red curve), still an asym- 
metry is observed. This originates from a rather broad 
distribution of the photoelectron energy. 

We can now analyze the dependence of the streaked 
lines upon various pulse parameters. Those with most 
influence on the streaking mechanism are the pondero- 
motive potential Up of the streaking field and the dura- 
tion of the pump pulse, tx- In Fig. |4] the dependence of 
the Auger line width of the Xe NOO transition is shown 
for a set of XUV pulse durations for a scaling parameter 
of 7 = 10. For larger pulse durations, a longer period of 
the slope of the streaking vector potential is accessible, 
consequently leading to a larger overall width, which is in 
accordance with the typical streaking mechanism. How- 
ever, the asymmetry with respect to the sign of the slope 
of the vector potential is more pronounced for shorter 
pulse durations (2.8 fs). This can be attributed to the 
fact, that for longer pulse durations, the line width is 
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FIG. 6. (color) Coincidence energy spectra for Auger and 
photoelectrons for selected time delays tx calculated by solv- 
ing Eqs. ([l| and S. The integrated photoelectron (Auger 
electron) distribution is plotted with red dashed (blue solid) 
lines. Parameters are the same as in Fig. [l] (Xe NOO with 
7 = 10). 



dominated by the streaking part and for shorter pulse 
durations the chirp becomes dominant, which will be dis- 
cussed in detail in Sec. |Vl 

For different ponderomotive potentials of the streaking 
field, shown in Fig. [5] a similar picture arises: the larger 
the ponderomotive potential, the larger is the streaking 
contribution leading to a relative decrease of the observed 
asymmetry. However, we note that for very small Up 
the line width asymmetry must vanish because of the 
vanishing vector potential. 



D. Auger electron and photoelectron coincidence 
spectra 

Additionally to the individual kinetic energy spectra 
of the Auger electrons, Eqs. (II]) and ^ also allow for 
the calculation of coincidence energy spectra of both in- 
volved electrons. This gives more detailed insight into 
their correlated motion. An example for the Xe NOO 
transition (7 = 10) is given in Fig. [6] for five selected 
time delays tx- AH spectra are dominated by a diago- 
nal line from top left to bottom right, which indicates 
an energy correlation between photoelectron and Auger 
electron. This is due to an energy exchange between both 
and governed by energy conservation. For the field-free 
cases (±144 fs) and the maximum of the vector potential 
(0 fs) a rather sharp spectrum is observed, whereas at 
the (approximate) zero transitions of the vector poten- 
tial (±60 fs) the streaking mechanism gives broad energy 
distributions, for both, the photoelectron and the Auger 
electron. 

In addition to the integrated Auger electron spectra 
(blue, solid lines), the corresponding photoelectron dis- 
tribution is plotted with (red) dashed lines. A careful 
inspection reveals, that the prominent asymmetry with 



respect to the slope of the vector potential (±60 fs) ob- 
served for the Auger electron, is not present in the pho- 
toelectron spectra. This indicates, that the photoelec- 
tron distribution carries no energetic chirp. A more de- 
tailed description and a simple picture for this are given 
in Sec. HVGl 



III. SEMI-CLASSICAL SIMULATIONS 

To proceed further and compare quantitatively with 
current experiments, it is crucial to take into account the 
3D geometry of the atom and the true time scales, i.e. to 
avoid the scaling procedure by 7. Since this is not possi- 
ble utilizing TDSE simulations, it is necessary to turn to 
a (semi-) classical description of FA-PCI including both 
electrons, the ion and the streaking field. Our classical 
method describing PCI is motivated by successful previ- 
ous models for the field- free case [52]. 

The classical dynamics of both electrons is governed 
by Newton's equations (nie = 1), 



rp(t) 
rAit) 



Fp{rp,rA,t) , 
FA{rp,rA,t) , 



(8) 



where the photoelectron (Auger electron) is denoted by 
index P (A). The propagation is split into two phases: (i) 
before Auger decay {tp. < t < tA^) and (ii) after Auger 
decay (i^^ <t <t^), where i^ is the time at the detector, 
which determines the corresponding forces in Eq. ([8|): 



• Phase (i) [tp^ <t <tA, 



Fp{t) = -VV+{rp)-EL{t) 



(9) 



Phase (ii) [t > t^J: 



Fp{t) =^ -yV^+{rp) - EL{t) -VVe-e{rA,rp), 
FA{t) = -VV^+{rA) ~ EL{t) -VVe-e{rp,rA).m 

For (i) only the photoelectron is propagated in the com- 
bined field of a singly charged ion, V^+(r), and the streak- 
ing field, EL{t). In phase (ii) both electrons experience 
the potential of a doubly charged ion, V'^~^{r), the streak- 
ing field, and their binary interaction Ve-e{i'i,r2). All 
interaction potentials are of pure Coulomb type: 



V+{r) 



1 



V^+ir) 



and Ve_e(ri,r2) 



■r-2| 



(12) 



By setting Ve-e = and additionally considering V^ = 
y^+, all e-e interactions and PCI effects can be turned off 
(denoted by "neglecting e-e interaction" in the following). 
The set of equations (l8| is completed by associated initial 
conditions 

'f'pitPi) ^Pp, and rp{tp^) = rp^ , (13) 

fA{tAi)=PA, and rA{tA,)^rAi ■ (14) 




electron at 
detector 



FIG. 7. Temporal parameters and electron distributions: the 
XUV pulse is centered at tx with a FWHM duration of tx- 
At ipj during the pulse the photoelectron is excited which 
triggers Auger decay at a time instant tAi ■ The measurement 
is performed long after the pulses and the decay are over, at 
time td (td — >■ 00). 



A. Initial condition sampling 

To reproduce the quantum mechanical nature of pho- 
toionization and Auger decay in our classical model, we 
developed a Monte Carlo (MC) sampling procedure for 
the initial conditions (131 and (14). During the XUV 
pulse, the photoelectron is released with the probability 



(proportional to the instantaneous intensity of the XUV 
pulse, ex E\) 



VrAn) 



1 



\J^TX 



exp 



(15) 



which creates the core hole at a time r^ = ip. — tx ■ The 
vacancy is filled after the time ta = tAi ~ ^Pi > by 
lifting the Auger electron into the continuum according 
to the decay law (probability density, see Fig. [t] for no- 
tations) 



VtaIta) = r^e 



-Tata 



(16) 



The kinetic energy distribution of the photoelectron fol- 
lows a Gaussian distribution. 



VeAEp.) 



1 



V2 



nax 



exp 



{Ep^ - Epg) 
2al 



(17) 



(11) with the spectral width ax centered around the energy 
Ep^ = (jjx — Ip, with the ionization potential of the core 



electron Ip. The (undisturbed) line shape of the Auger 



electron with mean energy Ea associated with Eq. ( 16 ) 
is a Lorentzian distribution 



VEAiEA,) 



Ta/27t 



lr2 

4^ A 



(18) 



(Ea, - Ea^ 
With that, the absolute values of the initial momenta are 



set by Eqs. (17) and (18) to 



IPpJ = y^2Ep: and |p^J = y^W^'^ 



(19) 



For small initial distances rp. and ta^ of the electrons 
from the ion, it is important to take into account the re- 
maining finite binding potential at the point of appear- 
ance of the electrons, V~^{rp.) and y^"'"(r^.), to assure 
their correct asymptotic momenta on the detector. En- 
tering as a free parameter in our model, we carefully 
checked the influence of different values of rp. and r^. 
ranging from 1 to 20 (in units of the Bohr radius) and 
found no significant change of the results. 

The directions of rp. and r^^. as well as of pp. and 
Pa. are given by the quantum mechanical angular dis- 
tributions of the associated initial state, approximated 

by 

^pmM = ^[1 + /3pm(3 cos V-l)] (20) 

with the asymmetry parameter /3, being available in the 
literature, e.g. 05]. We note as a technical aspect, that 
sphere point picking [35] is crucial for the correct MC 
sampling of Eq. ( 20 ) to maintain the correct uniform dis- 



tribution of points on a sphere. 



B. Extraction of observables 



We propagate Eq. (Sl) with initial conditions ([l3]) and 



(14 1, randomly distri 



Duted according to Eqs. 



utilizing a velocity Verlet algorithm with an adaptive 
time step size control, see e.g. [47]. This method will 
be called "Monte-Carlo Molecular Dynamics" (MC-MD) 
simulations in the following (MD refers to the classi- 
cal propagation of both interacting electrons leaving the 
atom). 

For each run, the final momenta pp and p^ of typ- 
ically 10^-10^ trajectories are recorded and sorted in 
angle- and energy-resolved histograms until convergence 
is reached. The Auger electron kinetic energy spectra 
are then obtained by integrating over a detector angle el- 
ement of 12.5°, typical for experiments, around the field 
polarization axis iz. Two opposite detection directions 
are possible, determined by the direction of A. We will 
only show results for the detector with positive energy 
shift at the maximum of the single-cycle vector poten- 
tial; the second detector gives the same results, but for 
changed sign in A. In experiments it is often favorable 
to consider two opposing detectors to assure the same 
streaking conditions [2QJ- Post-processing of the Auger 
line shapes is performed similar to the TDSE case, cf. 
Sec. |IIC| Additionally, as in the previous part, we restrict 
ourselves to the case of Auger electrons, the analysis of 
the photoclcctrons can be performed in a similar way. 



C. MC-MD-Results 

We may now drop the scaling procedure ([7| intro- 
duced for TDSE simulations and restore the true time 
constants. The result for a full scan of time delays for 
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FIG. 8. (color online) Auger line shapes of the Kr MNN 
transition for a set of time delays tx in a 1 THz single cycle 
streaking field with Up = 80meV for a XUV pulse duration 
of 28 fs at a photon energy of 97 eV. Results are obtained 
by MC averaging of MD trajectories. Shown is the case in- 
cluding photoelectron-Auger electron interactions (blue solid 
lines) and neglecting e-e interactions (red dashed lines). Note: 
the latter lines are shifted by 100 meV towards higher energy 
for clarity. 



the Krypton MNN transition in a 1 THz streaking field 
with a ponderomotive potential of 80meV is shown in 
Fig. Is] for the cases (i) including (blue solid lines) and (ii) 
neglecting (red dashed lines) e-e interactions. A similar 
picture as for the TDSE simulations, cf. Fig. [T] arises. A 
prominent asymmetry with respect to positive and neg- 
ative time delays, i. e. A{t)\t=tx — -^p > and Ap < 
respectively, can be found for (i) which completely van- 
ishes for (ii). Note: again, the lines with PCI effects 
excluded are shifted towards higher energy by 100 meV 
for better comparison. The shift is smaller compared to 
Fig. [T] due to the fact that 7 = 10 overestimates PCI in 
the case of the TDSE simulations. 

The FWHM and position of the line for the Xe NOO 
and Kr MNN decays are shown in Fig. [9] At the consid- 
ered photon energies of 91 cV for the former and 97 eV for 
the latter, photoelectron energies of 24 eV and 2.6 eV at 
comparable Auger electron energies of 34 eV and 40 eV 
are observed. Due to the slow photoelectron, for Kr a 
dramatic increase of PCI in comparison to Xe is expected, 
which is connected with a stronger chirp on the Auger 
electron's energy. This is confirmed by our calculations 
(red solid lines vs. blue dashed lines). If e-e interactions 
are neglected, similar line shapes and widths are observed 
for rising and falling fiank of the vector potential (black 
dotted line). These observations are in qualitative agree- 
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FIG. 9. (color online) FWHM of the Auger lines (right axis) 
of the Xe NOO and Kr MNN transitions calculated utiliz- 
ing MC-MD simulations. For Xe (Kr) the photon energy of 
91 eV (97 eV) leads to a photoelectron kinetic energy of 24 eV 
(2.6 eV). All other parameters are the sames as in Fig.^ The 
position of the maximum of the line, £max for both cases is 
shown in gray. The case neglecting e-e interactions for Kr is 
given by the black dotted line. 



ment with TDSE simulations discussed in Fig.lSland con- 
firm that PCI is the origin for the Auger electron's chirp. 
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FIG. 10. (color online) Line shapes of the Xe NOO transition 
(7 = 10) at three different time delays tx, corresponding to 
the FWHM maxima (see curve in Fig.lslfor 91 eV) and tx ~ 0, 
for streaking with Up = 98.8 meV. Results from TDSE simu- 
lations (top row) and MC averaged MD simulations (bottom 
row) are given. The cases including PCI (blue solid lines) 
and neglecting PCI (red dashed lines) are compared for falling 
(left), rising (center) and zero (right) slope of yl(i). 



D. Comparison with TDSE 

By construction, the MD simulations neglect any quan- 
tum effects in the electron dynamics, such as coherence, 
interference and spin. To test the above-introduced tech- 
nique, a detailed comparison of the line shapes calculated 
utilizing MD and TDSE methods for three different time 
delays is presented in Fig. Wo\ The Auger electron spec- 



trum of the Xe NOO transition, necessarily scaled by a 
factor of 7 = 10 for both simulations, is given for Ap < 
(left), Ap > (center) and Ap — (right) for situations 
including (solid lines) and neglecting (dashed lines) PCI. 
For better comparison, the Auger spectra obtained from 
TDSE and MD simulations have both been renormal- 
ized. This rescaling is necessary due to the small XUV 
ionization cross section, which has been neglected in the 
classical simulations. 

As a first observation, the line shapes obtained by MD 
simulations (bottom) are slightly broadened in compar- 
ison to the TDSE (top). This can be attributed to the 
averaging over the finite detector acceptance angle of 
12.5° in the 3D MD calculations. Here, trajectories are 
collected, which have been streaked with smaller ampli- 
tude due to their initial deviation (angular distribution) 
from the field-polarization axis. For both types of simu- 
lations, the line for Ap > is significantly broader than 



for Ap < which completely vanishes if PCI is turned 
off. Furthermore, both methods reproduce a similar PCI- 
induced shift of the line to higher energies. Thus the 
general trends as well as the underlying mechanism for 
the description of the asymmetry are correctly captured 
by the MD model and quantum effects in the electron 
propagation play no dominant role for the line width in 
the considered excitation regimes. 



IV. ANALYTICAL MODEL FOR AUGER LINE 
SHAPES 

In the previous sections, we have shown, utilizing 
TDSE and MD simulations, that Auger emission is 
chirped if PCI is involved which has a prominent impact 
on the line shapes in external laser fields. To get deeper 
insight in the underlying physics, we derive closed ex- 
pressions for the line shape of the Auger electron in the 
streaking field including PCI effects based on a classical 
ID model. 



A. Time-to-energy mapping 

The key mechanism of streaking is the mapping be- 
tween a temporal process and the measurable energy or 



momentum distribution. For Auger electrons, the tem- 



poral distribution follows the decay law, Eq. (16). The 
corresponding probability to find the Auger eli 
the continuum at a time t is given by 



ectron in 



VA{t)= / dTAVxArA) , 



(21) 



which approaches unity for long times (see Sec. Ill A and 
Fig. 7] for notations). The distribution (21) is translated 
by the streaking field to energy, thus the quantity of in- 
terest is the kinetic energy change of the Auger electron 
measured at a remote detector at time td- Its final mo- 
mentum is given by pAi + ^PitAi)- The field-induced 
momentum change evaluates to 



Ap : 



dt E{t) ^-A{tA,) 



(22) 



Ai 



where vanishing of the vector potential for t = t^ with 
td — >■ 00 is assumed. With that, we obtain for the Auger 
electron energy change 



E\ 



1 



kin 



[pA, ~ A{tAjf + AE 



PCI 



pX 



(23) 



A possible energy exchange between photoelectron and 
Auger electron due to post-collision interaction is ac- 
counted for by AE^^^{tx,TA)- It depends on the dis- 
tance from the ion, i. e. on the Auger time delay ta, and 
the pump-probe time delay tx- In the following we con- 
sider fixed (sharp) initial momenta of the two electrons, 
pp. and pAi ■ 

Let us first assume an infinitesimal duration of the 
pump pulse {tx — >■ 0), which corresponds to tp. = tx- 
An extension of the model to finite XUV pulse durations 



will be presented in Sec. IV F Expanding A around tp. 
to second order gives for the TA-dependent energy shift 



^PCl — 



Kn+PA.Ap: 



4'^^ ~ ~PA, ( ApTA 



Apri +Ai?^^Hr^). (24) 



Here, we use the notations Ap = A{tp.), Ap — 
d/dtA{t)\t=tp^ and Ap = d'^ /dt^A{t)\t=tp^ and neglect 
higher-order terms 0{t\uj\,A^). Eq. (24) translates 



the temporal distribution of Auger electrons governed 
by Eq. (21) to the energy domain through action of the 



streaking vector potential and PCI. This procedure was 
first applied in ,23j for the timc-to-energy transforma- 
tion due to PCI without external fields. In the present 
paper, we demonstrate, extending the simplified model 
of Ref. [20], this mapping including PCI and streaking, 
which gives direct access to closed expressions for the 
Auger line shape of FA-PCI. 



B. Lineshapes neglecting PCI 
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FIG. 11. (color online) Analytical model for the Auger line 
shape neglecting PCI, Eq. (26 1, compared to corresponding 
MD simulations (excluding PCI and sampling of pA^ ) at zero 
transitions of the vector potential for increasing (left) and 
decreasing (right) slope of A. Shown is Xe NOO in a 1 THz 
streaking field with a duration of 1 ps and a ponderomotive 
potential of 100 meV. The insets show the same data semi- 
logarithmically. 



probe time delays tx for zero transitions and maxima of 

A{t). 



1. Zero transitions of the vector potential. 

Since Ap ^ 0, the leading contribution to the map- 
ping (24) is linear in ta and higher-order contributions 



can be dropped, which gives 



TA 



PA.Ai 



(25) 



Substituting expression (25) in Eq. (21) gives the Auger 
line shapes for increasing (+) and decreasing (-) slope of 
A. 



/i±(e)=rie±^i% with Fi 



with the normalization conditions 



PAi\Ai 



(26) 



/O POO 

de/i+(e) = l, and / de/i_(e) = l. (27) 
-00 Jo 



The comparison with ID MD simulations for tx — > 0, 
neglecting PCI and without sampling of the initial mo- 
mentum pAi^ is shown in Fig. 11 for Xe NOO decay in 
a 1 THz streaking field. The streaked lines exhibit the 
same exponential decay law as the time dependence of 
the core hole decay. The direction of the slope of A only 
affects the orientation of the exponential tail. Deviations 
of Eq. ( 26 ) from the numerical solution are very small 



Let us first consider the case AE^'~'^ = in Eq. ( 24 ) and are due to the linearization of A and are only visible 



and find the Auger line shapes at characteristic pump- in the logarithmic representation (insets in Fig. 11) 
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FIG. 12. (color online) Auger lineshape for minima (left 
panel, Ap > 0) and maxima (right panel, Ap < 0) of A. 
Displayed is Eq. [29] compared to MD simulations. Parame- 
ters are the same as for Fig. |11| 



2. Extrema of the vector potential 

For maxima {Ap < 0) and minima [Ap > 0) of the 
vector potential, Ap — holds, thus the second order 
in TA is the leading contribution in Eq. (24). Because 



obviously ta > 0, we obtain only one solution 



2e 



TA 



ApPA 



(28) 



The corresponding line shapes for maxima ( "-" ) and min- 
ima ( "-f " ) oi A evaluate to 



/2±(e) 



Ta, 



PAi\Ap\ 



,(29) 



with the normalizations 

fO 



/U /"CJO 

de/2+(e) = l, and / de /2-(e) = 1 . (30) 

The comparison with MD data is given in Fig. [12] The 
line is dominated by a sharp onset at zero and a rather 
rapid decay. The second order expansion of A gives per- 
fect agreement with the simulation (logarithmic repre- 
sentation given in inset in Fig. 12). 



C. Analytical model for PCI 



We now consider the case AE^'^^ 7^ in Eq. (24) 



To obtain closed expressions for the lineshapes including 
streaking and PCI, the (semi-)classical model introduced 
in Sec. [Til] needs to be simplified in order to calculate 
^^pci Following [23l[24], we neglect the direct electron- 
electron interaction and model the PCI energy exchange 
by an instantaneous change in the ionic binding potential 
of X^+ to X+ for the Auger electron and X+ to X^+ for 




tp <ti< tp + r^ < t2 < t* < ts 



distance from ion r 

FIG. 13. (color online) Scheme of simplified ID propagation 
with PCI effects: in phase 1 {t < t*), the Auger electron (filled 
circles) catches up with the slow photoelectron (open circles) 
and overtakes it at i = t* (phase 2), pA > Pp is assumed. 
The propagation towards the detector in phase 3 (t > i*) is 
similar to phase 1. Figure after ref. |20| . 



the photoelectron. The propagation scheme, sketched in 



Fig. 13 then reads as follows: 



• Phase 1 {t < t*): propagation of the photoelec- 
tron (t > tp.) and the Auger electron (f > t^.) 
without interaction in the streaking field EL{t): 

tp^-El {t) and ta^-El (t) (31) 

with initial conditions 

rp{tp,) = rp., ppitpj = pp. 

rA{tA,) =rAi, PAitA,) ^PAi (32) 

• Phase 2 (t = t*): The Auger electron overtakes 
the photoelectron, changed screening of the ion's 
charge leads to energy exchange zbAii^^'-'^ — 1/r* 
corresponding to a momentum change of 

PP^Pp=PP + App ( A^pci , r ) 

PA^p\=PA + ApA{AE^^\t*) (33) 

• Phase 3 (t > t*): similar to phase 1 but with 
initial conditions 



rp{t*)^rA{t*)^r*; 
pp{t*) ^pp a.nd pA{t*) =Pa ■ 



(34) 



A straightforward integration of Eq. (31) gives the time 
of overtaking 



t* 



tA, 



TAPP - TAP 



(35) 



Pap 
and the corresponding distance from the ion 

r* = TAPr + Sr* . (36) 

Here, we introduced the notations 

pp =PP^ - M^P,) , 
PA =PA, - M^A,) , 
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FIG. 14. (color online) Field-assisted PCI without streaking 
contribution at zero transitions of A. Shown is Eq. (39 1 for 
A > (red dashed line) and A < (blue dotted line) in 
comparison to the field free case, Eq. ( |41[ ). Data is for N-00 
transition of Xe in a 1 THz streaking field with Up = 100 meV 
for Sr* = 1. 
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FIG. 15. (color online) Auger lineshapes at zero transition of 
A including streaking and PCI contributions. The analytical 
line shapes, Eqs. ( |45[ ) and ( |46[ l [bold, red lines] are compared 
to MD simulations (filled area) . The case of Ap > {Ap < 0) 
is shown in the left (right) panel. Parameters are the same as 
for Fig. [14] 



and 



rAP = rA, 



^P. 



AtA{i) , 



Pr 



PAP =PA-PP = PAi - PPi + MtPi) ~ M^Ai) 
PA.PPj 
PAP 

pATp, -ppTAj - A{tAi){rAi - rpj 

PAP 



With that, we obtain from Eq. (24) the ta dependence of 



the time-to-energy mapping function including PCI and 
streaking: 



e^ (ta) = -pA^ApTA 



1 



PrTA 



(37) 



The distance Sr* depends on the initial coordinates of the 
two electrons and their field-changed initial momenta. In 
most cases 5r* will be a small correction to r* . However, 
for situations with slow photoeletrons and relatively fast 
Auger electrons, i. e. situations with strong PCI, Sr* 
may become large. In the following, we derive general- 
ized Auger line shapes for FA-PCI, improving the results 
presented in [20|, where Sr* = was assumed. 



D. FA-PCI without streaking 

Before using the full mapping, let us neglect the first 
term in Eq. (37) linear in ta that is attributed to the 



streaking contribution discussed before. Then we have a 
hyperbolic mapping function 



,PC1 



1 



PrTA 



(38) 



Utilizing Eq. ( 38 1 , the straightforward transformation of 



the time distribution (21) gives for the PCI-induced en- 



ergy change of the Auger line shape 

f3±ie) = r,\e^^^'^'~'/^\ withFg^^ (39) 
e-^ Pr 

f.l/i5r* 

de/3(e)-l, (40) 



where "-I-" ("-") refers to ip > {Ap < 0) inpr and Sr* . 
From Eq. ( 39 1 we can immediately read off the energy 



distribution for the field-free case (and assuming Sr* 
0), 



/f(e) 



111 



-rA/ip^e) 



(41) 



with pjl = pPiPAi/{pPi ~ PAi)i in accordance with the 
result given in [23 . Our result ( p9| ) differs in the way, that 
although we exclude the explicit streaking contribution 
in Eq. (37), the field-changed initial momenta pp and 
Pa are included. An example of Xe NOO in a 1 THz 
streaking field is shown in Fig. 14 In addition to the 



case of positive and negative slope of A, the field-free 
case, Eq. (|4T|), is displayed. For this specific case, no 



strong influence of the field on the pure PCI process is 
visible. However, f^^ has a slightly higher maximum, 
corresponding to smaller width, in contrast to the effect 
observed in the simulations in the previous section (note: 
this asymmetry is not the observed chirp), /g is exactly 
in the middle between both. 
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E. Auger line shapes including FA-PCI 



Using the full mapping function ( 37 ) gives a quadratic 
equation for ta, 



ta+'^a 



Pre 



PCI 



5r*pAiAi 



6r* 



^PCI 



-1 



PrPA.Ap 



PrPA,Al 



= . (42) 



For the inversion of Eq. (37), we assume Pri^A) = Pn 



hence we neglect the additional implicit ta dependence, 
which enters through the vector potential. We define 



e± = e± —or ; a = ^; 

4 2paAp 



P ^ ''J^^^ (43) 



and obtain 



TA± = — ae+ ± ak^ 



(44) 



with k± — \ (r_ ± |/?|. For Ap > 0, only the positive 

branch ta^ can be realized {ta > 0), which gives for the 
lineshape with r4 — lalF^: 

k: — e 



/+(£) = r4 



-r4(fc+-e+) 



(45) 
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FIG. 16. (color online) Time-to-energy mapping of the Auger 
electrons in the streaking field for Ap > (left) and Ap < 
(right). The mappings with out PCI, Eq. 1 25 1 blue dashed 



lines, without streaking, Eq. (38 1 green dotted lines, and in- 
cluding PCI and streaking, Eq. ( 37 1 red solid lines, are shown. 



Difi'erent signs of Ap lead to a drastic change of the mapping: 
for Ap > the whole energy space is accessible in a bijective 
way, whereas in contrast for Ap < a forbidden region below 
£0 occurs and "early" and "late" Auger electrons are mapped 
to the same energy. 



For Ap < 0, both solutions (44), ta+ and ta~, are 



possible. Thus, the temporal distribution function (21 ) is 



split into two parts, /g^""" e ^^^^dTA + J^. e ^-^'^-^dr^, 
where Ti„i„ separates both branches, ta+ and ta^, at 
e(Tinin) = cq. The straightforward transformation of both 
integrals to energy gives for the joint energy distribution 
function 



/_(e) = 2F4e 



-r4e+ 



coshF4fc_ — sinhF4A:_ 



(46) 
The line shapes of Xe NOO for both cases, Ap > and 



Ap < 0, are given in Fig. 15 for 6r* = 1. While for the 



former, the line is broadened by PCI, for the latter, the 
line is compressed and completely different line shapes 
for subsequent zero transitions of A with different sign of 
the slope are observed. As in the previous cases, perfect 
agreement with simulations based on numerical solutions 
of Eqs. (3Tp4| by means of MC averaged MD simulations 



(in analogy to Sec. Ill but without momentum averag- 
ing) is observed, and the linearization of A has, in the 
considered regimes of pulse duration and ponderomotive 
potential, no significant influence on the streaked Auger 
spectra. 

To explain the strikingly different shape of the Auger 
lines in Fig. |15[ the mapping functions from time to en- 
ergy are shown in Fig. [16] for the same set of parameters. 
The linear streaking part contributing to Eq. ( 37 ) is plot- 



ted with blue dashed lines, the hyperbolic PCI term with 
green dotted lines and the sum of both results in the red 
solid lines. By comparing Ap > (left) and Ap < 
(right), the cause for the different lineshapes becomes 



visible: whereas for the former, both terms add up to 
a bijective mapping function spanning the whole energy 
axis from — oo to oo, for the latter one a forbidden en- 



ergy region for e < eg = 2\ pA\Ap\/pr occurs (gray line 
in Fig. 16 1. This leads to a drastic compression of the line 
(right panel in Fig. 15 1, where Auger electrons released 



at two different time moments can be mapped into the 
same energy interval. This situation is completely ab- 
sent for Ap > which leads to a broad distribution of 
Auger electrons, cf. 



15 This effect is 



left panel in Fig. 
a direct consequence of the interplay between the hyper- 
bolic PCI-induced chirp on the Auger electron energy, 
AE^'^^ ex 1/r* oc 1/ta, and the linear "chirp" intro- 
duced by the streaking field, where the sign of the latter 
depends on the direction of the streaking field at the time 
of the core hole creation. An experimental verification of 
this mechanism utilizing XUV pulses from FLASH and 
HHG exciting the Xe NOO and Kr MNN transitions has 
been presented in Ref. |20j . The comparison of the ex- 
perimentally obtained Auger electron spectra with the 
theoretical results calculated based on MC-MD simula- 
tions, as presented in Sec. Ill shows perfect agreement. 



F. Finite XUV pulse duration 



Eqs. (45) and (46 1 describe the shape of the Auger 



energy distribution for infinitesimal pulse duration of the 
XUV excitation. To account for finite pulse durations 
Tx , a similar transformation of probability distributions 
from time to energy as for the case of the decay law 
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FIG. 17. (color online) Auger lineshapes for zero transitions 
of A neglecting PCI contributions for a finite XUV pulse dura- 
tion of 20 fs FWHM. All parameters refer to THz streaking of 
the Xe NOO transition. The convoluted line shape, Eq. (51 1 
[bold red line] , is plotted against MD simulations (blue area) . 
The case of infinitesimal XUV pulse duration, Eq. (26 1, is 
shown also for comparison (blue dashed lines). 
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FIG. 18. (color online) The same as Fig. 17 but for the case in 



eluding PCI. The result of a numerical convolution, Eq. (52 1, 
red bold lines, is compared to MD simulations with sharp 
momenta (blue area) and with proper sampling over initial 
momenta pp. and pAi (blue dashed lines). Additionally, 3D 
MD simulations, cf. Sec. Ill are shown (orange dotted lines). 



(211 needs to be performed. The temporal distribution 



of a photoelectron excited by a Gaussian shaped pulse is 
described by Eq. (15 1. At zero transitions of the vector 



potential, utilizing the same linearization of A as was 
used for Eq. ([24|, we obtain for the Tj-dependent energy 
shift (see Fig.TTffor notations) 



l{n) = -p,A{tx)n + 0{A^,Tf) 



(47) 



Using this mapping function, a streaked energy spectrum 
due to the finite XUV pulse duration can be calculated: 



fx{i) 



1 



hrTx- 



■ exp 



^ XI 



(48) 



with the normalization condition 

difxii) - 1, and Txi = PAAitx)Tx ■ (49) 

At each photoelectron "birth" time tp., the Auger 
"clock" starts, and with that the energy mapping of the 
temporal distribution of Auger electrons. Therefore, the 
final streaked Auger energy distribution is given by the 
convolution 



/oo 
difx{e)fA{e-i) 
-OO 



(50) 



If the PCI contribution is neglected, i. e. 
cf. Eq. (26), the integration in Eq. (50 1 



fA{e)^h±{e), 
can be carried 



out analytically and gives for the line shape 



fxi±ie) 



2 



ir?e±i^^'erfc (^^^^ 






± 



XI 



(51) 

The lineshape for two subsequent zero transitions of A 
for an XUV pulse duration of 20 fs FWHM is shown in 
Fig. 17 (bold red line). The finite excitation time interval 



of the core hole broadens the pure Auger decay line (blue 
dashed line) . Despite the rather long pulse duration com- 
pared to the core hole lifetime, the exponential decay of 
the case of infinitesimal excitation duration is still im- 
printed on the convoluted line. Corresponding MD data 
(blue area, according to method (i) below) is in perfect 
agreement with Eq. (51). 

Considering the PCTdistorted line shapes, /yi(e) = 
/±(e), the integral 



fx±{e) 



defxi~e)f±{e-e) 



(52) 



cannot be solved analytically. The result of a numerical 
integration of Eq. ^ for a 20 fs FWHM XUV pulse is 



given in Fig. 18 for positive and negative slope of A{t). 
Although the asymmetry with respect to Ap > and 
Ap < is less pronounced than for the case of infinites- 
imal XUV excitation duration shown in Fig. [T5j still a 
difference between ascending and descending slope of A 
is visible indicating the chirp in Auger emission. To verify 
the accuracy of the analytical result, additionally three 
different sets of simulations are shown: (i) numerical so- 
lutions according to Eqs. ( 31p4 ) [(blue) filled area], (ii) 
similar to (i) but with proper averaging over the initial 
momenta pA^ and Pp^ [blue dashed lines] and (iii) 3D MD 
simulations according to the scheme in Sec. |III[ also in- 
cluding angular distributions [orange dotted lines], (i) 
resembles the assumptions of the analytical model, ex- 
cept for the linearization of A, and shows perfect agree- 
ment with Eq. (52). Solutions according to (ii) and (iii) 



show a substantial broadening of the Auger line. For (ii) 
this results from the natural Auger line width and the 
bandwidth of the XUV pulse, and for (iii) in addition 
from the angle integration. This broadening occurs in 
a similar way for rising and falling fiank of A and does 
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FIG. 19. (color online) Energy (left axis) and FWHM (right 
axis) of the Auger electron (dotted) and the photoelectron 
distribution (PCI: solid, without PCI: dashed). Shown is the 
case of Xe NOO scaled by 7 = 10 obtained within TDSE 
simulations. Parameters are the same as for Fig. [4] 



not affect the observed asymmetry attributed to the PCI- 
induced chirp. Therefore, the analytical line shape model 
(52 1 catches the important features of FA-PCI. It can be 



evaluated numerically for a large set of parameters due to 
its simple convolution structure. Thus, Eq. ( 52 ) is well- 



suited for the detailed investigation of the properties of 
FA-PCI and its dependence on the streaking conditions 
and XUV parameters. 



G. Photoelectron distributions 

In the previous sections, we identified the physical 
mechanism for the observed time-dependent chirp on the 
Auger electron's energy: a direct connection between the 
time instant of decay and an associated (unique) energy 
shift mediated through post-collision interaction. At this 
point, a remark on the consequences for the correspond- 
ing photoelectron distribution is appropriate. As a mat- 
ter of fact, the kinetic energy of the photoelectron is af- 
fected in a similar way as it is for the Auger electron, 
but, by reason of energy conservation, with an opposite 
sign. Thus, the photoelectron is slowed down due to PCI 
by the same amount of energy ^E^''^^ the Auger electron 
has gained. 

However, this does not result, as one might guess, in a 
chirp on the photoelectron energy distribution with dif- 
ferent sign, as already pointed out in Sec. |IID| The re- 
sults of TDSE simulations, carried out as described in 
Sec. [Ill are shown in Fig. [19| The FWHM for the photo- 
electron line is depicted for the case including PCI (solid 
line) and neglecting PCI (dashed line) for a full set of time 
delays between pump and probe pulse. For both cases, no 
asymmetry with respect to the rising and the falling flank 
of the vector potential is observed. Only a broadening of 
the line for the PCTincluded case is present, resulting in 
an equidistant upward shift of the PCI curve in compari- 
son to the case neglecting PCI, and therefore, no chirp on 
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FIG. 20. (color online) Auger line shape of the Kr MNN 
transition in a 3.3 THz streaking field at fixed ponderomotive 
potential of 100 meV for different XUV pulse durations with a 
photon energy of 97 eV. Shown are solutions of the analytical 
model, Eq. ([52 1. 



the photoelectron energy can be identified. The promi- 
nent asymmetry in the FWHM for the Auger electron is 
shown for comparison with dotted lines in Fig. |19[ 



Returning to Eq. ( 37 ) , the mapping from the time mo- 



ment of Auger decay, r^, to energy, the origin for the 
absence of a chirp becomes clear: while for Auger elec- 
trons, each decay time corresponds to a certain amount 
of energy transfer, for photoelectrons no such connection 
can be found. For each time moment of photoemission, 
every Auger decay time ta is possible, and with that any 
arbitrary energy transfer due to PCI. Thus, only a PCI 
broadening of the line is expected, which is the same for 
every release time of the photoelectron and explains the 
observed delay dependence of the FWHM of the photo- 
electron spectrum in Fig. [19) 



V. USE OF FA-PCI FOR PULSE 
CHARACTERIZATION 

In this last part we ask the question, whether the PCI- 
induced asymmetry with respect to the direction of the 
slope of the vector potential may help to improve the 
pulse characterization capabilities of the streak camera. 
In Sec. [m we already commented on the influence of the 
ponderomotive potential and the XUV pump pulse dura- 
tion on the observed asymmetry. This raises the question 
whether this dependence, especially on the pulse duration 
Tx, can be used as a sensitive parameter to recover the 
pulse durations in experiments. Due to the larger effect 
of PCI, we choose Kr MNN decay at a photon energy of 
97 eV in the following. 

The line shape for different XUV pulse durations at 
the zero transitions of the vector potential with Up = 
100 meV is shown in Fig. [20| calculated according to 



Eq. (52). The left panel shows the broadened line / 
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for ^ > and the right panel the corresponding com- 
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FIG. 21. (color online) The same as Fig. 20 but for different 
ponderomotive potential of the streaking field at fixed XUV 
pulse duration of 28 fs. 



pressed line f^ for A < 0. For the shortest pulse dura- 
tions (3.5 and 7fs) with tx ^ T^ the largest asymmetry 
is observed, whereas for long pulses (42 fs) no clear dis- 
tinction between A > and A < is possible. The 
largest impact on this asymmetry has the compressed 
line, which is sharp in the case of very short pump pulses, 
cf. Fig. 15 right panel. For increasing tx, the convolu- 



tion with the Gaussian shaped time distribution of the 
XUV excitation smoothens (broadens) this line, until the 
streaking-induced broadening predominates. This result 
agrees qualitatively with the TDSE simulations presented 
above in Sec.|TrCl cf. Fig. [I] 

A similar behavior is observed upon change of the pon- 
deromotive potential of the streaking pulse at a fixed 
pump pulse duration of 28 fs, see Fig. [21] Here, the 
largest asymmetry is observed for small Up, whereas the 
asymmetry gradually decreases upon increase of Up. This 
effect can be as well explained by a domination of the 
streaking contribution for large Up, where the broaden- 
ing of the line due to the larger momentum transfer from 
the streaking field exceeds the PCI contribution. Again, 
a similar picture arises in TDSE calculations, cf. Fig. [5) 

In a further step, we want to evaluate the asymmetry 
in more detail. To this end, we introduce a classification 
parameter ^, defined as 



\w — w^\ 



(53) 



where w^ is the FWHM of the Auger line at A > and 
A < 0, respectively. This parameter ^ G [0, 1] describes 
the relative asymmetry and is zero for vanishing asym- 
metry and attains finite positive values smaller than one 
in any other case. 
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FIG. 22. (color online) Position of the Auger line's maximum 
Em, the corresp ond ing width w and the asymmetry param- 
eter ^, cf. Eq. (53 1, vs. XUV pump pulse duration. Data is 



for Kr MNN decay at a ponderomotive potential of 100 meV 
of the streaking pulse. The forbidden region for Ap < at 
infinitesimal short XUV pulse duration is marked by eo- The 
superscript "-{-"("-") refers to yip > (Ap < 0). 



A. Effect of XUV pump pulse duration 

A key property for successful pulse characterization is 
the existence of observables sensitive to the XUV pulse 
duration. In Fig. [22] several properties of the Auger line 
shapes are plotted depending on the XUV pulse dura- 
tion at fixed Up — 100 meV: the position of the maxi- 
mum of the line, e^, the FWHM of the line, w^, and 
the corresponding asymmetry parameter ^. Let us first 
consider the energetic position of the line. For A < 0, 
e^ is positive, stemming directly from the positive tail of 
the decay law. For small pulse durations, the forbidden 
region e < eo, discussed in Sec. |IVE] is observed with the 
sharp onset of the line at eq. With increasing pulse dura- 
tion, the position of the maximum shifts towards higher 
energies, a clear consequence of the convolution with the 
Gaussian shaped temporal distribution of the XUV pulse. 
This shift saturates for large pulse durations, tx S> T^ , 
due to the broad convolution function. For A > and 
e^, a similar trend is observed for e+ , but with decreas- 
ing energy of the maximum. Due to the PCI-induced 
broadening of the line, the effect is here less pronounced 
than for the compressed line. 

The width of the line, w^, shows the typical strong 
infiuence on the XUV pump pulse duration: the larger 
Tx, the larger w^. This phenomenon is the basic prin- 
ciple of the streak camera utilized for the estimation of 
pulse lengths. Here, the observed asymmetry due to the 
PGl-induced chirp on the Auger electron's energy mani- 
fests itself in two branches for the width: one for A > 
labeled by w~^ and one for A < indicated by w^ , cf. 
dashed lines in Fig. [22] For conventional chirp-free situa- 
tions both branches coincide (no asymmetry with respect 



16 



0.5 



-? 0.25 

3 



-0.25 



- 










£,« 


^^^,„„— -^^r^^'*'^ 






\ 


-—■ --'' 


_,'-'' 




^^^^""^ 


..■-■<w^ .. 


,.'''''W~ 




^ 






\ ^ 


^ ,,.--''' 


,--''' 






^y''' ,..-- 




- 


/'/ 


1 







1.5 
1.25 
1 
0.75 



0.5 X 



0.25 



50 100 

ponderomotive potential Up [meV] 



150 



FIG. 23. (color online) The same as Fig. 22 but for varied 



ponderomotive potential Up of the streaking field at a fixed 
pump pulse duration of 28 fs. 



to i > and i < 0). 

Since with two opposing detectors, both situations 
^ > and A < can be recorded simultaneously, also 
the asymmetry parameter ^ carries valuable information 
about the single-shot pulse properties. Its pulse-duration 
dependence is plotted by the (green) dotted line: for 
short pulses a large asymmetry of about 0.75 is observed 
with a rapid decrease down to about 0.1 at pulse dura- 
tions of 20 fs. The largest variation is found for pulses 
below 10 fs where tx is comparable to the Auger decay 
time r^ . Thus from measuring the width of the Auger 
lines at opposite slopes of the streaking vector potential 
simultaneously, a reconstruction of the pulse duration is 
possible, even if PCI effects are present. 



B. Ponderomotive potential of streaking pulse 

An important question is the dependence of the asym- 
metry on the streaking conditions and, in particular, the 
ponderomotive potential Up of the streaking field. This 
parameter is, in principle, easily tunable in experiments, 
either through the frequency ojl (limited by th e pu lse 

the 



23 



duration tx <C i^^^) or the intensity £■[. In Fig. 
position of the maximum, the FWHM and the asymme- 
try parameter for a scan of C/p, by variation of E^, at a 
fixed XUV pulse duration of 28 fs are shown. For vanish- 
ing Up, no streaking occurs, which gives vanishing energy 
shifts ei; and widths w^. [Note: no natural line widths 



are included in Eq. (52)] 



For increasing Up, the position of the maximum of the 
line shifts towards higher energies for the case of ^ < 
(e~) and towards smaller energies for 4 > (e+). This 
effect is similar to the behavior observed upon variation 
of Tx- Also, as in the previous case, two branches of the 
width can be identified, originating from the asymmetry 
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FIG. 24. (color online) Asymmetry parameter ^ for different 
ponderomotive potentials of the streaking pulse as a function 
of the XUV pulse duration. Parameters are the same as for 
Fig. [221 



transfer from the probing field to the electron increases 
and, with that, the streaking, resulting in broader lines. 
Again for the case of chirp-free XUV excitation without 
PCI effects the "-I-" and the "-" branches would coincide. 

Additionally, the asymmetry parameter ^ is plotted in 
Fig. |23[ Starting from a rather high value at very small 
Up, it exhibits a rapid drop when Up increases to about 
lOmeV, followed by a slow convergence with only a weak 
dependence on Up over a wide range of approximately 
100 meV. 

Fig. [24] summarizes the possibilities for pulse charac- 
terization using FA-PCI. The asymmetry parameter ^ is 
shown as a function of tx for different Up of the streak- 
ing pulse. For all considered values of Up a monotonic 
behavior with large asymmetry for small pulse durations 
and vice versa is observed. The larger Up the faster is 
the drop of ^ at small pulse durations and, with the cor- 
responding strong variation with tx, a high sensitivity 
of ^ in the range of pulse durations below 10 to 20 fs 
occurs. Thus increasing Up allows to extend the region 
of sensitivity to slightly larger XUV pulse durations. In 
conclusion, measuring ^ at a given ponderomotive po- 
tential of the streaking field for ^ > and A < simul- 
taneously utilizing opposite detectors, an estimation of 
the XUV pulse duration based on (timc-rcsolved) Auger 
electron spectroscopy is possible. The highest sensitivity 
is reached for XUV pulse durations below 10 fs with a 
rather strong variation of the measured parameter ^ by 
a factor of four. 



VI. CONCLUSIONS AND OUTLOOK 

In this paper, we gave a detailed theoretical explana- 
tion of the experimental observations in |20j , which show 
evidence of an energetic chirp in Auger emission. Based 
on solutions of the TDSE we could reproduce the chirp 
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for model systems and explain its origin by post-collision 
interaction. This formed the basis for classical model- 
ing of the photoelectron and the Auger electron in the 
continuum, including all electron- electron and electron- 
ion interactions. Using Monte-Carlo averaged Molecular 
Dynamics simulations for the electrons, we were able to 
verify this chirp. The quantitative comparison with ex- 
periments including detector resolutions and acceptance 
geometries using the Xe NOO and the Kr MNN transi- 
tions presented in [20] shows perfect agreement between 
our approach and the light-field driven streak camera in 
the considered range of parameters. For deeper insight 
and to obtain a more flexible tool, we derived a classical, 
analytical line shape model for the Auger electron that 
fully includes the XUV pulse shape, streaking and PCI 
effects and thus captures all important properties. We 
further showed, how our results may be used as a tool 
for estimating the length of XUV pump pulses if PCI 
effects are involved. 

In the present work, we focused on the Auger electrons, 
which was motivated by currently available experiments. 
The corresponding photoelectron distribution was briefly 
discussed, and we explained why no energetic chirp is 
present there. A detailed analysis will be part of a future 



work. Worthwhile considerations include the influence of 
additional XUV pulse parameters such as chirp and sub- 
structures, e. g. spikes as present in the case of free elec- 
tron laser sources. Additionally, it will be advantageous 
to extend our purely classical model for the FA-PCI line 
shape to account for quantum effects in order to describe 
interference and spin effects. 

Finally, it would be very interesting to investigate in 
experiments with either Kr MNN or Xe NOO the be- 
havior when the XUV photon energy is increased. If the 
proposed FA-PCI scenario is correct, then the chirp of 
the Auger spectra should vanish when the photoelectron 
energy starts to exceed the Auger electron energy. 
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